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The following changes should be noted: 


Figure 17, Eber’s equation should be: 


Nu = 0.0149 Re 0,8 

instead of 


Nu = 0.0217 Re°/ 8 


The line in figure 17 indicated as result- 
ing from Eber's equation should be changed 
to correspond to the change in the equa- 
tion. 
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HEAT— TRANSFER AND BOUNDARY-LAYER TRANSITION ON 

A HEATED 20° CONE AT A MACH NUMBER OF 1.53 

By Richard Schemer, William R. Wimbrcrw, and 
Forrest E. Gowen 


SUMMARY 


Heat— transfer data from supersonic wind-tunnel tests of a heated 
20 ° cone have been compared with theoretical results obtained by two 
methods for determining the convective heat transfer In laminar 
boundary layers in a compressible fluid. The cone was heated elec- 
trically and was tested at a Mach number of 1.53* Local rate of heat 
transfer and surface— temperature measurements were made over a range 
of Reynolds numbers and nominal surface temperatures with both laminar 
and turbulent boundary layers. 

The theoretical and experimental results in the case of the 
laminar boundary layer were found to be in good agreement in terms 
of the heat— transfer coefficients In the region on the test body 
where the theory was considered applicable. Good agreement In terms 
of rate of heat transfer was obtained by the use of the theoretical 
heat— transfer coefficients and the true temperature potential. 

The effect of heat transfer on boundary-layer stability wae indicated 
by surface-temperature measurements for a uniform power input distri- 
bution, the sudden decrease in surface temperature at the beginning 
of the turbulent boundary— layer region being indicative of the 
transition. The results provided a qualitative verification of the 
effect of heat transfer on laminar boundary— layer stability that 
had been predicted theoretically by Lees. (NACA Technical Note 
No. 1360.) 

The general heat— transfer equations developed In NACA TN 
No. I 3 OO are shown to reduce, for cones, to simple relationships, 
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and these are presented in the form of design charts by which the 
local rate of heat transfer may he determined on cones with 
attached how waves. 


INTRODUCTION 

Because of aerodynamic heating, the practical operation of 
aircraft at high speeds is dependent on the provision of adequate 
insulation and cooling systems for the aircraft structure, equipment, 
pay load, and occupants. The design of such systems. In turn, is 
dependent on the existence of adequate heat— transfer data end on 
the development of theories hy which the data may he correlated and 
Its application extended. 


The most extensive experimental Investigation to date In the 
field of heat transfer at high velocities was conducted in Germany 
hy Eher. (See reference 1.) This work provides the basis for most 
heat— transfer calculations for proposed supersonic aircraft. How- 
ever, the air flow In the test section of the supersonic wind tunnel 
at Kochel, in which Eher performed his experiments, was such that 
there has been some question as to the extent of the laminar boundary 
layer on the test bodies. (See fig. 5 of reference 1.) Since there 
are large differences in the rates of heat transfer through laminar 
and turbulent boundary layers, additional experiments have been 
needed to clarify Eber’s results. 


Another aspect of the heat— transfer problem, both at subsonic 
and supersonic speeds, is the effect of heat transfer on the stabil- 
ity of a laminar boundary layer. The theoretical work of Lees 
(reference 2) indicates that the effect of surface heating is 
destabilizing to a laminar boundary layer and also Indicates that 
the effect of surface cooling is stabilizing. The results presented 
in reference 3 for a very low Mach number are In agreement with 
the results of reference 2j however, no experimental data are 
available to indicate the effect of heat transfer on boundary-layer 
stability at supersonic speeds. 


The purpose of the Investigation presented In this report was to 
obtain heat-transfer data on a body of revolution with first, a laminar 
boundary layer and, then, a turbulent boundary layer, and to compare 
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these data vith the theoretical results calculated by the methods 
of references 4 and 5 and with the results obtained by Eber. The 
qualitative effect of heat transfer on the stability of the laminar 
boundary was also to be determined. 


SYMBOLS 


The following symbols have been used in the presentation of the 
theoretical and experimental data: 

A area, square feet 

a speed of sound, feet per second 

Cp specific heat at constant pressure, Btu per pound, °F 

cy specific heat at constant volume, Btu per pound, °F 

c arbitrary constant 

g gravitational constant, 32.2 feet per second squared 
H total pressure, pounds per square foot, absolute 
h local heat— transfer coefficient, Btu per hour, square foot, °F 
K average heat— transfer coefficient, Btu per hour, square foot, °F 

k thermal conductivity, Btu per hour, square foot, °F per foot 

l body length, feet 

M Mach number, dimensionless 

m Mach number parameter (~ My 2 ), dimensionless 


Nu average Nusselt number 
Nu local Nusselt number 


hs'N 
*s/' 

(S> 


dimensionless 


dimensionless 


Nug boundary— layer Nusselt number 


h&\ 

kJ’ 


dimensionless 


Pr Prandtl number x 3600g^, dimensionless 


p static pressure, pounds per square foot, absolute 
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total rate of heat transfer, Btu per hour 

local rate of heat transfer, Btu per hour, square foot 

gas constant for air, 1718 foot squared per second squared, °F 

/p v Ys\ 

Reynolds number l — J, dimensionless 
radius of body, feet 

distance from nose along surface of the body, feet 
tempera ture,°F absolute 

recovery surface temperature, °F absolute 

pseudo-surface temperature [ T a ’ = p (T Q — T v ) + T v ] , °F absolute 

fluid velocity parallel to the surface at any point within the 
boundary layer, feet per second 

fluid velocity Just outside the boundary layer, feet per second 
distance normal to the body surface, feet 

surface— temperature parameter; for a Prandtl number of 1 . 0 , 

/T '-T \ 

0 = ( £ Y J and for a Prandtl number of O. 73 , p = 

\T q — (T 7 / 

dimensionless 

ratio of specific heats (cp/c v ), dimensionless 

boundary— layer thickness, feet 
cone half -angle, degrees 

absolute viscosity, pound— second per square foot 
air density, slugs per cubic foot 
air density ratio (p/p a ), dimensionless 
unit surface shear, pounds per square foot 


SdEi) 
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In addition, the following subscripts have been used: 
a reference air density 

s fluid conditions at the body surface 

v any point along the body. Just outside the boundary layer 

x location of a particular limit of integration along the length 

of the body 

o fluid conditions at total temperature and pressure (after isenr- 
tropic compression from static conditions) 


1 fluid conditions Just behind an attached oblique shock wave from 
the nose of a body 

The superscript 1 together with the subscript s have been 
used to indicate the pseudo-surface temperature, T g *, and the 
physical constants of air based on this temperature, p. g ' and k g ’. 


ANALYSIS 

In order to obtain continuity in this report, the various theo- 
retical developments involved in the presentation and explanation of 
the test data are presented separately in appendices. Only the results 
of each development are presented in the text. 

A method for calculating the rate of heat transfer in the laminar 
boundary-layer region of bodies of revolution in steady supersonic 
flight is presented in reference 4 and is used as the basis of the 
theoretical calculations for the present investigation. The method 
assumes a linear velocity profile within the laminar boundary layer 
and also assumes a Prandtl number of one, but considers the effect of 
compressibility. The general equations of reference 4 are shown in 
appendix A to reduce, for cones, to the single equation, 


Nu = rr = 1 * 225 Re (A15) 1 

The equation defining the variable B, as a function of Mach number 

The equation designation (A 15 ) indicates equation (15) of appendix A. 
This method of designation is used throughout this report. 
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and surface— temperature parameter, is given in appendix A. 

Equation (A15) gives the value of local Nusselt number or local 
heat— transfer coefficient at any point on a cone. However, in many 
heat— transfer problems the average heat— transfer coefficient is 
required rather than the local value. The average value of heat- 
transfer coefficient from the nose to any point on a cone with a 
laminar boundary layer is shown in appendix B to be given by the 
relation. 


h = J h (B5) 

This simple relationship, first recognized by Hantzsche and Wendt 
(reference 5 ) results from the form of equation (A15) and the geometry 
of cones. This equation may be used to convert local values of 
Nusselt number to average values as long as the surface temperature 
is constant. 


In a laminar boundary layer in subsonic flow, the velocity pro- 
file is known to be very similar to the profile calculated by Blasius. 
Velocity profiles in a laminar boundary layer in supersonic air flow 
have not been measured for any appreciable range of Mach numbers, but 
the profiles have been calculated by several investigators. The trend 
of the calculated profile shapes with increasing Mach number is from 
the Blasius profile at subsonic Mach numbers toward an almost linear 
profile at a Mach number of 10. (See reference 6 . ) The effect of 
surface cooling at any Mach number is to make the velocity profile 
approach that of some lower Mach number, or to become less linear. 
Although a linear velocity profile is assumed in the development of 
the method of reference 4, the effect of this assumption is shown 
by the comparison between the methods of references 4 and 7 devel- 
oped in appendix C. The method of reference 7 assumes a Blasius 
velocity profile in an incompressible fluid and assumes a Prandtl 
number of one. Since the only differences in the two methods are 
the profile assumptions and the consideration of compressibility in 
the method of reference 4, the difference in the results obtained 
by the two methods at some subsonic Mach number, at which compressi- 
bility can be neglected, would only be due to the velocity profile 
assumptions. In the comparison of the heat— transfer coefficients on 
a flat plate, given by the two methods, the method of reference 7 
gives the relation 


h 


= 0.332 


and from the method of reference 4, for M=0 and 0=1.0, 


(C5) 
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h 


0.286 



(C7) 


It is evident from the constants in the above equations that the 
effect of the linear velocity profile assumption, at the conditions 
of zero Mach number and zero heat transfer, is to decrease the cal- 
culated heat— transfer coefficient by about 15 percent relative to 
that obtained by the method of reference 7« The assumption of a 
linear velocity profile leads to a calculated boundary— layer thick- 
ness that is 50 percent greater than is the case with the method of 
reference 7 and a boundary— layer Nusselt number that is 30 percent 
greater. These two effects are partially compensating and the 
difference in the heat— transfer coefficients, as indicated by the 
constants in equations (C5) and (C7), is relatively small. 

The local values of Nusselt number on cones with attached bow 
waves can be calculated by equation (A15) i.f the conditions of the 
air stream just outside the boundary layer are known. The details 
of the method by which the theoretical data based on reference 4 
were calculated for this report are presented in appendix D together 
with a step-by-step outline of the method for using a series of 
design charts based on equation (A15) . 


APPARATUS AND TEST PROCEDURES 
Wind Tunnel 

The tests were performed in the Ames 1— by 3 — Toot supersonic 
wind tunnel No.l. This tunnel was temporarily equipped with a 
1— by 2— l/2— foot test section and a fixed nozzle that provided a 
test— section Mach number of 1.53* Since no aerodynamic forces were 
to be measured, the strain— gage balance equipment was removed and 
the test cone was mounted with a suitable adapter to the balance 
housing. 


Test Cone 

The usual case of heat transfer at supersonic speeds is for 
heat to flow into the surface rather than out of the surface. From 
the theoretical aspect either case would be satisfactory to obtain a 
partial check on the theory of reference 4, but a complete comparison 
requires the testing of both a heated and a cooled body under similar 
test conditions. An electrically heated cone was chosen for these 
tests because of the simplicity of the experimental techniques which 
could be employed. 
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The 20 ° cone was constructed as shown in figure 1. The exterior 
shell was machined from stainless steel and all other metal parts 
were made of copper# The exterior surface of the model had a smooth, 
ground finish, estimated to be a 30— microinch root mean square (rms) 
surface. The walls of the shell were tapered to maintain an approxi- 
mately constant incremental resistance along the cone length when cold. 
The cone was heated by passing a high amperage (800 amperes maximum), 
low voltage (0.45 volts maximum), alternating electrical current 
longitudinally through the cone surface. Because no current would 
flow through the extreme nose of the cone, the forward 25 percent of 
the cone was in effect unheated. 

Eight thermocouples were installed at equal length increments 
along the cone to allow determination of the temperature distribution. 
The thermocouples were made from 30-gage copper-const ant an duplex wire 
with welded junctions. They were installed in holes drilled completely 
through the shell and were soldered in place. Ten leads of 20— gage 
copper wire were also installed in the shell, in a similar manner, to 
provide a means of measuring incremental voltage drops along the cone. 
The locations of the thermocouples and the voltage leads are indicated 
in figure 1. A photograph of the assembled cone is shown in figure 2, 
and a photograph of the cone installed in the wind tunnel is shown in 
figure 3. 


Instrumentation 

The wiring of the test cone was connected as shown in figure 4. 

The variable voltage transformer controlled the input to the primary 
side of the power transformer. The secondary side of the power trans- 
former was grounded to the tunnel shell which acted as one lead in the 
circuit. The other lead consisted of two parallel cables that were 
connected to two binding posts at the base of the cone. Two cables 
were used to keep the cable size down to a convenient diameter. These 
cables passed through a current transformer which was in turn connected 
to an ammeter to measure the current input to the cone. 

The eight thermocouples were connected through a selector switch 
to a potentiometer. The potentiometer was used to obtain a zero read- 
ing on an external light— beam galvanometer, the potentiometer output 
then being equal to the thermocouple potential. 

The ten voltage leads from the cone were connected through a 
selector switch to an electronic voltmeter in such a manner as to 
measure the voltage drops of successive increments along the cone. 

The local power input, or rate of heat transfer per unit length, 
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is given by the product of the current and the incremental voltage 
drop. 


The total temperature of the air stream was measured by nine 
thermocouples in the tunnel settling chamber which were connected 
through a selector switch to a direct reading potentiometer. 


Procedure 

Data were obtained over a range of Reynolds numbers from approxi- 
mately 0.5 to 2.5 millions. This variation of Reynolds number was 
effected by varying the pressure level within the tunnel. The tunnel 
was first brought to the desired pressure and then allowed to run 
until the general temperature distribution on the cone came to equilib- 
rium. When this condition was reached, the surface temperature of the 
cone was measured by the surface thermocouples. The surface temper- 
ature measured under these conditions (zero heat flow) is called the 
recovery surface temperature, or just recovery temperature Tr. 

The heating circuit was then closed and the cone heated to the desired 
temperature, as indicated by the potentiometer reading of the most 
forward surface thermocouple, by adjusting the input voltage. Since 
the average total temperature of the air stream was in the order of 
100° F, cone temperatures of 120°, l4o°, l60°, l80°, and 200° F were 
arbitrarily chosen as nominal values at which to obtain data. The 
surface temperature varied along the length of the cone through a 
range of about 5° bo 35°> depending on the temperature level, the front 
of the heated section of the cone always being the hottest. 

With the cone at the desired temperature, the following data 
were read and recorded: the total pressure and total temperature of 

the air stream, the current input to the cone, the incremental voltage 
drops, and the local surface temperatures of the cone. These data 
were obtained at each of the nominal cone temperatures previously 
mentioned and at nominal values of total pressure of 3, 6, 9 , 15 , 
and 21 pounds per square inch absolute. 

Upon completion of the tests described, surface roughness was 
employed to obtain data with a completely turbulent boundary layer. 
Approximately the first 2 inches of the nose of the cone were sprayed 
with clear lacquer and, before the lacquer was completely dry, it was 
sprayed again with lampblack in suspension in lacquer thinner. After 
the thinner evaporated, the lampblack adhered to the lacquer base and 
provided a band of fairly uniform roughness around the nose of the 
cone. Liquid— film tests were performed to determine if the roughness 
was sufficient to cause premature transition. It was found that at 
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total pressures above 6 pounds per square inch absolute, the boundary- 
layer was completely turbulent. Tests similar to those previously 
described were performed at nominal total pressures of 9 and 15 pounds 
per square inch absolute with the completely turbulent boundary layer. 


ACCURACY OF RESULTS 


The accuracy of the experimental data was determined by esti- 
mating the uncertainty of the individual measurements which entered 
into the determination of the final results. The over— all uncer- 
tainty of any given parameter was then obtained by geometric summa- 
tion of the uncertainty of each of the factors entering the final 
value of that parameter as indicated by the method employed in 
reference 8. 

The estimated uncertainty of the basic measurements are as 
follows : 

Total temperature T 0 ±2° F 

Recovery surface temperature Tp ±0.5° ? 

Free— stream temperature just outside the boundary layer . T v ±2° F 
Surface temperature Ts ±0.5° F 


Total pressure H Q ±0.05 centimeters of mercury 

Incremental voltage drops ^ -2 percent 

Input amperage I±8 amperes (l to 3 percent) 

Cone dimensions inch 

Cone segment surface areas .±3*2 percent 

The calculated accuracy of the final parameters are as follows: 

Surface— temperature parameter P±2.^ percent 

Temperature potential AT ±4.3 percent at P=1.4 

±1.5 percent at (3=2.0 

Local rate of heat transfer, q ±4.1 to ±5.0 percent 
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Nusselt number 


Heat— transfer coefficient 


.h ±4.4 to ±6.6 percent 
Nu +4.4 to +6.6 percent 


Reynolds number 


Ee ±1.8 to ±1.9 percent 


A further error was introduced in the experimental data by radi- 
ation of heat energy from the cone to the tunnel walls and an effort 
was made to determine the order of magnitude of the radiation by 
experimental means. The total heat transferred from the cone is 
equal to the sum of the convective and radiant heat transfer plus the 
end losses, and the radiant heat transfer is proportional to the 
difference in the fourth powers of the absolute temperatures of the 
cone anri wind— tunnel wall. The convective heat transfer is a function 
of total pressure and will become zero when the total pressure is 
reduced to zero. Therefore, at zero total pressure, the heat transfer 
will be entirely due to radiation. Since it is impossible to evacuate 
the tunnel to zero pressure and measure the heat transferred by 
radiation directly, the heat loss due to radiation was evaluated from 
the data obtained at the various test conditions with the tunnel in 
operation. 2 

The total heat transferred Q as measured at the various 
pressures was divided by the difference in the fourth powers of the 
cone and tunnel-wall absolute temperatures, and the resulting parameter 


pressures. Because the surface temperatures along the cone were not 
equal, data obtained by cross— plotting was used in the determination 
of the correction for radiant heat transfer. The ordinate of the 
logarithmic plot at zero pressure is a measure of the heat transferred 
to the tunnel walls by radiation and includes such factors as the 
Stephanr-Boltzman constant, the shape factor, and the emissivities of 
the cone and walls. The quantity thus attained, however, was so 
small as to be completely masked in the ±5~percent uncertainty of 
the measured heat transfer. Solutions for several of the elements 
gave slightly negative losses. Consequently, the correction for 
radiation was assumed to be negligible. 

2 An attempt was made to obtain the radiation calibration with the 
tunnel inoperative, but the cone surface temperatures were found 
to be very erratic because of free— convection currents. For this 
reason the method was abandoned. 

3 The method of reducing the data to constant values of surface- 

temperature parameter is discussed in detail in the section of this 
report titled "Results and Discussion." 



was plotted logarithmically against the corresponding total 
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Conduction along the akin of the cone also affected the data at 
the hase and at the nose. Calculations show that about 10 percent 
of the total heat generated in the first heated element is conducted 
to the unheated nose portion and a slightly higher percentage is lost 
from the last element through the base of the cone. Data from the 
first and last elements have been neglected in the analysis of the 
test data and the elements between these two appear to receive as 
much heat from neighboring elements as they lose. Consequently, the 
conduction losses are assumed to be negligible. The test Mach number, 
1.53, was selected as the average of the linear Mach number gradient 
in the region in which the model was installed and the maximum devi- 
ation from the average Mach number was approximately +0.02. The Mach 
number gradient in the test section was neglected in the reduction of 
the test data. 


RESULTS AND DISCUSSION 

The measurements of local power input were converted to local 
rates of heat transfer by dividing by the incremental areas and 
converting the electrical units to heat units. Heat— transfer coeffi- 
cients were obtained from the local rates of heat transfer by dividing 
by the temperature potential (T a — Tp). Nusselt numbers were obtained 
by the combination of the appropriate values of heat— transfer coeffi- 
cient, reference length, and thermal conductivity as previously 
defined. 


Laminar Boundary— Layer Heat Transfer 

The surface— temperature distributions along the cone for various 
nominal values of surface temperature are shown in figure 5* The 
temperature variation with length is due to the local values of 
electrical resistance and the heat— transfer— coefficient distribution. 
The heat— transfer rates in a turbulent boundary layer are, in general, 
much greater than those in a laminar boundary layer; therefore, the 
sudden drop in surface temperature toward the base of the cone, which 
appears in figures 5(<i) and 5(©), is indicative of transition to 
turbulent flow within the boundary layer. 

A rigorous comparison of theory and experiment would require 
constant values of surface— temperature parameter and hence constant 
surface temperature along the length of the cone. However, the surface 
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temperatures obtained in the experiments were not constant. In order 
to make the desired comparison, the original data were plotted as local 
heat—' transf er rate as a function surface— temperature parameter, as 
shown in figure 6. Lines for each longitudinal segment of the cone 
were drawn through the data points. The values of local heat-transfer 
rate for constant values of surface— temperature parameter were then 
obtained by cross— plotting. 

The comparison of the data, on the basis of constant values of 
surface— temperature parameter with a changing surface temperature, 
involves the assumption that the small variation in surface temper- 
ature ahead of a particular point on the cone does not affect the 
heat transfer at that point. The validity of this assumption is 
illustrated in figure 7 by a comparison of data for a nominal surface 
temperature of l80° F and cross— plotted data for a surface— temperature 
parameter of 1.8 on the basis of the Nusselt number — Reynolds number 
relationship. The lines in figure 7 are faired through the data 
obtained on the aft portion of the cone only. The difference between 
the two methods of data presentation is small, and for this reason 
the comparison of theory and experiment in terms of constant values 
of surface— temperature parameter is valid for the present experiments. 

It should be noted that the cross— plotted data are indicated by flagged 
symbols. This method of indicating cross— plotted data has been used 
throughout this report. 

The effect of the large variation in surface temperature which 
occurs at the beginning of the heated portion of the cone (s/z=0.25) 
is indicated by the initially decreasing values of Nusselt number 
with increasing Reynolds number for each tunnel pressure. (See fig. 

7.) This effect can be explained by consideration of the changes which 
occur in the boundary— layer temperature profile as the layer flows along 
the cone. 

The local rate of heat transfer at any point on the cone is 
given by the product of the thermal conductivity of the air adjacent 
to the surface and the slope of the boundary— layer-temperature profile 
at the surface. At the beginning of the heated portion of the test 
cone, a relatively cold boundary layer flows onto the heated area 
and the slope of the boundary— layer-temperature profile becomes large 
because of the large difference between the air and surface temperatures. 
The air temperature at the surface will approach the surface temperature 
as the air continues to flow along the heated surface, or the local rate 
of heat transfer downstream of the surface— temperature discontinuity 
will approach the value that would have existed if the surface- 
temperature discontinuity had not been present. The data Bhown in figure 7 
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indicate this tendency toward the rear of the cone where the change 
in surface temperature with length is small in comparison with that 
near the nose . 

The effect of heating the aft portion of the cone will he to 
increase the boundary— layer thickness in this region. The resulting 
laminar boundary— layer thickness can be calculated, at least approxi- 
mately, by the following method: If the surface temperature is 

assumed to be discontinuous at the edge of the heated region (no 
longitudinal conduction), the boundary— layer thickness at any point 
on the cone, either heated or unheated over its entire length, can 
be calculated by the method of reference 4. The boundary— layer 
thickness at any point along the heated portion can be approximated 
(as shown in fig. 8) as the thickness of the boundary layer for a 
completely heated cone less the difference in boundary— layer thicknesses 
at the edge of the heated portion, for a completely heated cone and 
for an unheated cone. The correction obtained by this method is small 
at the beginning of the heated region on the cone, and, because the 
correction is small, it can be expected to be reasonably accurate 
at any downstream position. At some point far downstream, where the 
boundary— layer thickness is considerably greater than at the beginning 
of the heated region, the percent error in boundary— layer thickness 
would be insignificant. 

A comparison of the theoretical local heatr-transfer coefficients 
for the heated— and unheated— nose conditions with the experimental 
values for a surf ace— - temperature parameter of 1.4 is shown in 
figure 9. The agreement between the experimental and the theoreti- 
cal values, corrected for the effect of the unheated nose, is good 
over the after portion of the cone. Figure 9 also indicates the 
failure of any method for calculating heat— transfer coefficients, 
based on boundary— layer thickness, when a large change occurs in 
the assumed relation between the boundary— layer velocity profile 
and temperature profile. A method based on different assumptions 
is needed to calculate the local rates of heat transfer in regions 
where large surface— temperature gradients exist. Such a method will 
be necessary in order to calculate the optimum location of surface- 
cooling heat exchangers for high-speed aircraft. 

The experimental and theoretical values of local Nusselt number 
are shown as functions of length Reynolds number and surface- 
temperature parameter in figure 10. The theoretical values are 
corrected for the effect of the unheated nose by the method illustrated 
in figure 8. The effect of the correction is to alter the slope of 
the lines from the 2:1 slope indicated by equation (A15) . The 
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correction also produces the discontinuities between the lines for 
each value of total pressure, because the value of boundary— layer 
thickness to which the Nusselt number is related is not directly 
related to the length Reynolds number for the different values of 
total pressure. 

It should be noted in figure 10 that the agreement between slopes 
of the theoretical lines (reference 4) and the trend of the data is 
almost exact within the small scatter of the data. This agreement is 
indicative of the accuracy of correcting the boundary— layer thickness 
for the effect of the unheated nose, because the thickness correction 
primarily affects the exponent of the Nusselt number to Reynolds 
number relationship. 

The relationship derived by Hantsche and Wendt (reference 5) 
for laminar boundary layers is also plotted in figure 10 for comparison 
with the experimental data and the comparable results calculated by 
the method of reference 4. 

For a rigorous comparison of theory and experiment, the effect 
of each of the test conditions on the final results should be known. 

The following variables affect laminar boundary-layer thickness and 
therefore the heat transfer at any point on a test body: 

1. Distance along the body 

2. Velocity of air flow along the body 

3. Ambient— air temperature 

h. Surface temperature 

5. Ambient— air pressure 

6. Surface— pressure gradient 

7. Surface— temperature gradient 

8. Surface roughness 

The effects of the first five of these variables are accounted 
for in equation (A15), and the effect of surface-pressure gradient 
has been eliminated from the experiments by the selection of a cone 
for a test body. However, a small pressure or Mach number gradient 
does exist in the wind-tunnel nozzle. An approximate correction 
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for the effect of surface-; temperature gradient on boundary— layer 
thickness has been made in the comparison of theory with experiment. 
However, the effect of the surface-temperature gradient on the 
assumed boundary-layer temperature and velocity- prof ile relationship 
cannot be included in the correction. The effect of surface rough- 
ness on laminar boundary— layer thickness is not known quantitatively 
for the surface finish on the test cone. Since surface roughness 
will be present to some extent on all supersonic aircraft, its effects 
should be investigated, at least to the extent of determining a value 
of roughness below which there will be little or no effect on laminar- 
boundary— layer thickness or stability. 

The agreement between the theoretical results based on references 
4 and 5 and the experimental data, shown in figure 10, is satisfactory 
over the rear portion of the cone where the theories are considered 
to be applicable. Whether or not the comparison is favorably or 
adversely affected by surface roughness, pressure gradients or surface- 
temperature gradients can only be determined by fundamental investiga- 
tions of each of these effects. 

The plots of experimental local heat-transfer coefficient against 
surface-temperature parameters from which cross plots were made indicated 
a slight decrease in heai>-transfer coefficient with reduction in surface- 
temperature parameter; however, this trend was within the range of the 
experimental accuracy (±6 percent) as is the trend indicated by the 
theory of reference 4 (±k percent). For this reason, the heat- 

transfer-coefficient distributions of figure 11 are shown as only 
functions of total pressure. 

Satisfactory agreement between theory and experiment is shown 
in figure 11 over the rear portion of the cone where the theory is 
considered to be applicable. Because of this agreement, in terms of 
heat— transfer coefficient, poor agreement between the theory of refer- 
ence 4 and experiment in terms of local rate of heat transfer can be 
expected because of the incorrect temperature potential in the theo- 
retical equation [q=h(T s *— T 0 ) ] . Therefore, it appears logical to use 
the true temperature potential, corresponding to a Prandtl number of 
0.73 (T s -Tr) in the theoretical calculations of local rate of heat 
transfer rather than that corresponding to a Prandtl number of one 
(Ts'-Tq) that a rigorous interpretation of the theory would dictate. 

The desirability of this empirical change in the theory is indicated 
by the more satisfactory agreement between theory and experiment 
shown by the curves for a Prandtl number of 0.73 in figure 12. 
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Laminar-Boundary— Layer Stability 

The streamwise extent of the laminar boundary layer on the 
cone for various surface temperatures is indicative of the effect 
of heat transfer on boundary— layer stability. The distance from 
the nose of the test body to the transition point was obtained from 
the surface— temperature distribution curves of figure 5. The inflec- 
tion points on the curves of figure as indicated by the black dots, 
were selected as being the average transition points. The effect of 
heat transfer on boundary— layer stability is shown in figure 13 as a 
plot of the length Reynolds number at the transition point against 
the average surface— temperature parameter up to the indicated tran- 
sition point. The decrease in transition Reynolds number with 
increasing surface— temperature parameter confirms the prediction of 
reference 2 and agrees with the experimental results of reference 3« 
The curve of figure 13 shows a hyperbola— like variation of the tran- 
sition Reynolds number with surfaces temperature parameter, indicating 

/'cKeA 

that the rate of change of boundary— layer stability J decreases 

with decreasing stability. The difference between the values of 
Reynolds number for transition at 15 and 21 pounds per square inch 
total pressure is believed to be due to a change in air— stream turbu- 
lence level. 


Turbulent Boundary-Layer Heat Transfer 

The experimental surface temperature and local rate of heat- 
transfer distributions along the 20° cone with an artificially 
induced turbulent boundary layer are shown in figures lb and 15 . 

The same data in the nondimens ional form of local Nusselt number 
as a function of length Reynolds number are shown in figure 1 6 . 
Because of the scatter, a line of 0.8 slope (the known slope for 
turbulent boundary layers) was faired through the data points. 

The points which are displaced farthest above the line are those 
from the forward portion of the cone and, as in the case of the 
laminar boundary— layer data, are affected by the surface— temperature 
discontinuity at the beginning of the heated region. 

A comparison of the average values of Nusselt number from the 
turbulent boundary— layer data with the results obtained by Eber 
(reference l), by Hantzsche and Wendt (reference 5)> and those 
obtained for a laminar boundary layer from the design charts 
(appendix D) and corrected by the four— thirds factor to obtain 
average values of heat— transfer coefficient (appendix B) is made 
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in figure 17. The curves for each nominal value of surface temper- 
ature in figure 17 tend to approach a common line asymptotically. 

This result occurs because the percent effect of the large local 
Nusselt numbers in the nose region (due to the surface— temperature 
discontinuity) gradually decreases as more of the cone is included 
in the average. The asymptotes of the experimental lines were drawn 
with a slope of 0.8 which is also the slope of the line given by 
Eber's equation. 

It is evident from a comparison of the various curves of figure 
17 that Eber’s results were obtained from test bodies with turbulent 
boundary layers. The low Beynolds numbers of Eber’s tests should 
have produced laminar boundary layers; therefore, the transition must 
have been caused by external disturbances. Figure 5 of reference 1 
shows the great number of shock waves which existed in the test section 
of the Kochel supersonic wind tunnel in which Eber conducted his 
experiments. It is known that such shock waves are very effective 
in causing premature transition of the laminar boundary layer. The 
fact that transition was induced artificially in both Eber’s and the 
present experiments limits the applicability of the data. The dif- 
ference between the turbulent boundary— layer data from the present 
experiments and the results given by Eber's equation is probably due 
to the difference in the methods of causing transition. In the present 
experiments transition was induced by roughness at the nose of the cone 
and the boundary layer was entirely turbulent . In Eber's experiments 
transition, due to shock waves, would be expected to occur farther aft 
on the cone and the boundary layer at the nose would be laminar. 

This being the case, the average heat— transfer coefficient and the 
average Nusselt number obtained by Eber should be lower than those 
obtained in the present experiments. The scatter of the data obtained 
by Eber would have masked any change in the slope of the Nusselt 
number — Beynolds number line that would be expected to result from 
mixed laminar and turbulent flow. 

It follows from the preceding discussion that any turbulent 
boundary— layer heat -transfer data which are not obtained with natu- 
ral transition or knowledge of the preceding laminar boundary layer 
will not be generally applicable to the calculation of the cooling 
requirements of supersonic aircraft. The fact that Eber’s equation 
gives usable results when applied to the specific problem of calcu- 
lating the temperature— time relationship of the skin at the nose of 
missiles indicates that turbulent boundary layers exist in this region 
or that the method of calculation rather than the data determines the 
results obtained. 

If turbulent boundary layers do exist in the nose region of 
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missiles, improvements in shape and surface condition would allow 
longer runs of laminar boundary layer with the result that the rate 
of increase of surface temperature with time would he materially 
reduced. Also, the heat capacity or mass of the skin could be 
reduced for a given rate of increase of surface temperature. This 
latter effect would provide an improvement in mass ratio, and, there- 
fore, an improvement in the range of the missile. The weight advan- 
tage of maintaining laminar boundary layers to reduce the required 
capacity of aircraft cooling systems would be apparent with any 
method of cooling. 


CONCLUSIONS 

The results of the foregoing investigation lead to the follow- 
ing conclusions: 

1. Experimental heat— transfer coefficients obtained from tests 
of a heated 20° cone at a Mach number of 1.53 have been found to be 
in satisfactory agreement with two theoretical methods of calculating 
the rate of heat transfer in the laminar boundary— layer region of 
bodies of revolution in a compressible fluid. 

2. Satisfactory agreement was obtained between the theoretical 
rates of heat transfer based on NACA TN No. 1300 and those determined 
by experiment, in the region of the test body where the theory is 
considered applicable, when the theoretical heat— transfer coefficients 
and the true temperature potential were employed. 

3. The theoretical prediction of Lees (NACA TN No. 1360) that 
the effect of heating a surface with a laminar boundary layer to a 
temperature above the recovery surface temperature is to destabilize 
the boundary layer, has been confirmed experimentally at a Mach number 
of 1.53. 


Ames Aeronautical Laboratory, 

National Advisory Committee for Aeronautics, 
Moffett Field, Calif. 
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APPENDIX A 

LOCAL HEAT-TRANSFER COEFFICIENTS ON CONES 

The general equation for laminar boundary— layer thickness on 
a body of revolution in a compressible fluid is, from reference 
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It should he noted that the physical properties of the air in the 
preceding theoretical equations and in the following equations for 
a Prandtl number of one are referred to the psuedo-surface temper- 
ature T s 1 . This change in the nomenclature from reference 4 is 
necessary for the comparison of theory and experiment on the basis 
of equal values of surface— temperature parameter but for different 
values of Prandtl number. 

For the more specific case of a cone, the surface— pres sure 
coefficient is constant for a given Mach number, and equation 
(Al) becomes, 
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(A2) 


Because the radius of a cone is a linear function of its length, the 
integral of equation (A2) reduces to a constant ( 1 / 3 ), times the 
length ratio — , or, at a given point on the surface 

u 


3 l 



(A3) 


The Reynolds number for the flov just 
layer using the viscosity based on surface 


outside the boundary 
temperature is 


Ee 


Pv 7 s 
V-s 


(aU) 


Also, since in the experimental investigation the value of air- 
density ratio a* vill be one, it can be eliminated from equation 
(A3). With these simplifications, the laminar boundary— layer thick- 
ness relation for any cone becomes 


-2 _ 2 s2 
” 3 BEe 

or 

0.8l6 s 

o = — — 

V^BEe 


(A5) 
(A 6) 


The expression for surface shear per unit area for the linear veloc- 
ity profile of reference k is 
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Reynolds analogy between skin friction and heat transfer for 
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compress ib le flow gives the relation 

la = -Q _ 

V c p (T s -T E ) 
q , 

V 

c pM- 

Since Pr = -g— = 1.0, by assumption, then = T 0 and 

= ^ g »(T^ -T 9 ) 

* v 


(A8) 


(A9) 


(A10) 


Combining equation (A10) with equation (7) gives the relation 

kg' . 

<J = 4- (Ts'-T 0 ) 


and since 


then 


q = h (T s '-T 0 ) 


Wt = 1.0 = Nu & 

K s 


(All) 


(A12) 


(A13) 


With the laminar boundary— layer— thickness relation and the boundary- 
layer Nusselt number relation known, the two can be combined to give 
values of local heat-transfer coefficients directly. 


V/Blte 

0.8l6s 


or in terms of local Wusselt number 


(Aik) 


Nu = ka_ - 1.225 Vb Re 

k' 


(A15) 
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APPENDIX B 

AVERAGE HEAT— TRANSFER COEFFICIENTS 

The local rate of heat transfer in the laminar boundary— layer 
region of a cone can be expressed in equation (Al4) as 




(Bl) 


and the incremental area over which the local heat-transfer coeffi- 
cient is applied can he shown to he 


dA = c 2 s ds 


(B2) 


The average value of the heat— transfer coefficient is then given by 
the relation R 


or 


/ h dA ci (cb) //s ds 
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but, since 


then 


h = Cij/Js 


h = ^ h 


(B5) 


This relation has also heen obtained, in a slightly different form, 
by Hant&sche and Wendt in reference 5* 
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APPENDIX C 


COMPARISON WITH INCOMPRESSIBLE-FLOW THEORY 

For a flat plate, the laminar boundary— layer— thickness relation 
of reference 7 reduces to the form 


6 


2 5«3s 2 

Re 


(Cl) 


for the boundary— layer thickness measured at the point in the 
velocity profile where the dynamic pressure is one— half of free— 
stream dynamic pressure. Also, in reference 7, it is shown that for 
the Blasius velocity profile the boundary— layer Nusselt number is 
given by the relation 


Nu 5 = ijE. = 0.765 


(C2) 


The following relation is obtained from the method of reference 4 
for the boundary— layer thickness measured at the same point in the 
velocity profile: 


5 2 = 2s ‘ 


B Re 


(C3) 


and 


Nu 6 = — = 1.0 
& k ' 


(C4) 


Combining and rearranging equations (Cl) and (C2) gives the relation 


h = 0.322k /- 1 

S^ 


(C5) 


Similarly, equations (C3) and (C4) yield the relation 


h = 




(C6) 
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For a Mach number of zero, and zero heat transfer (0 = 1.0), k = kg* 
and 

h = 0.286k | (C7) 

A comparison of the constants of equations (C5) and (C7) is indicative 
of the effect of the linear— velocity— profile assumption. 


APPENDIX D 

METHOD OF CALCULATION AND DESIGN CHARTS 

The value of the surface— temperature parameter can be calculated 
from the known boundary— layer conditions by the relation (for Prandtl 
number = 0.73)* 


0 


0.73 


VjJV 

V*v 


(Dl) 


where 

Tr = T v (l+Jto 2^ V^) (D2) 

With a Prandtl number of unity as is assumed in the theory 
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1-0 


V-^v 

T o“*v 


(D3) 


In order to have similar temperature profiles in the actual and 
theoretical cases, the surface— temperature parameters must bo equal. 

Olo-Oo.t* W 

Therefore , the pseudo-surface temperature is given by the relation 

T s ' = 3(T 0 ^V) + T v (D5) 


or 
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T s '=?^ + T v (D6) 

because 

J* - i]Sr (D7) 

x o V 

With the values of surface— temperature parameter and pseudo-surface 
temperature known, the values of the parameter B, the viscosity, 
and the thermal conductivity of air at the surface based on the 
pseudo-surface temperature can be determined. This, in turn, allows 
the Reynolds number corresponding to the desired position on the 
cone to be calculated 


Re 


PyVS 

Ms' 


(D8) 


With the values of B and Reynolds number known, the local Nusselt 
number can easily be determined by equation (A15). The local heat- 
transfer coefficient can be determined from the local Nusselt number 
by the relation 


h 


= Nu 


V 


(D9) 


The theoretical results presented were calculated from the 
foregoing relations. The conditions of the air stream just outside 
the boundary layer were obtained by the use of reference 9, rather 
than by the linearized theory of reference 10 as indicated in 
reference k. With this change, the limit of applicability of the 
method is not the extreme body fineness ratio dictated by linearized 
theory, but rather the Mach number for nose shock-^wave detachment. 

The change in limiting fineness ratio requires the length s in the 
foregoing equations to be taken as the slant length because the 
assumption in reference 4 that the surface and axial lengths are 
equal is not valid for blunt bodies. 


The following outline gives a step— by— step procedure for 
calculating the rate of heat transmission to a cone moving at 
constant supersonic velocity. Use is made of the charts of this 


CONFIDENTIAL 


NACA EM No. A8L28 


CONFIDENTIAL 


27 


report (fig. 18) which, were developed from the theory set forth in 
reference L. Table II of reference 9 is very useful in many of the 
calculations, and its application is indicated in the appropriate 
steps. However, the symbols used in reference 9 differ from those 
used in reference 4 and the present report. A table of equivalent 
symbols follows: 


Present report 
and 

reference 4 

Reference 9 

T/T 0 

T/T a 

Ty/To 

T/T a at M = My 

Pv/ Po 

p/p a at M = My 

a v/ a o 

a/a a at M = My 


To begin the calculations the following information must be 
known: 

M flight Mach number 

T ambient static air temperature, °F absolute 

o 

T s surface temperature to be maintained on the cone, F absolute 
6 C half-angle of the cone, degrees 
p ambient-air pressure, pounds per square foot 

The calculations then proceed with the determination of the following 
parameters : 

1. Total temperature T 0 



or, enter table II of reference 9 with M and find T/T 0 directly. 
( See the preceding table for equivalent symbols . ) 
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2. Mach number Just outside the boundary layer of the cone My.. 
Enter figure 18(a) with M and 0 C and determine My/M. (This figure 
is taken from reference 11.) 

3« Temperature of the air stream j.ust outside the cone boundary 
layer T v , 



or, enter table II of reference 9 with M = My and find T v /T 0 
directly. 

4. Recovery surface temperature Tp. Enter figure 18(b) with 
My and determine Tr/T q . 

5. Surface— temperature parameter 3, 


0 = ± 302 . 

Tr-tv 

6. Pseudo-surface temperature T g ', 

T s « = 3(T 0 -T y ) + T v 

7. Total pressure behind bow shock wave Hi. The total pressure 
ahead of the bow shock wave H Q is given by 



or, enter table II of reference 9 with M and find p/H directly. 
Then 

7 l 


Hi 

(7+ l)M 2 sin 2 9 

7-1 

2y M 2 sin 2 0— (7— l) 

Ho 

-(y-l)M 2 sin 2 0+2^ 


L 7+1 -■ 


where 9 is the bow shock— wave angle and can be determined from 
figure 7 of reference 9« 

8. Density at total pressure behind the bow shock wave p 0 . 
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= Hi_ 

P oi ET 0 


9- Density just outside the cone boundary layer p v 


Py 


P 


°1 



or, enter table II of reference 9 with M = My and find p y /p 0i 
directly. 

10. Velocity of sound at total temperature conditions a Q 

a 0 2 = ?KT 0 


11 . 

the cone 


Velocity of air stream just outside the boundary layer of 

V, 



or, enter table II of reference 9 with M = M^ and find 
V = My X a v . 


a v/ a o 


then 


12. Absolute viscosity at the surface of the cone 
figure 18(c) with (T s ’— 460) and determine p s *. 


|i s ’. Enter 


13. Beynolds number per foot of slant length Ee/s. 


Be = P^V 
s p s » 


14. Eeynolds number for various positions on the surface of 
the cone Ee. 

(a) Choose stations along the surface of the cone at 

which it is desired to determine local heat— transfer 
rates. 

(b) Measure the distances s along the surface of the 

cone from the apex to the stations in feet. 
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(c) Then, the Reynolds number is equal to 
each station. 



x s for 


15. Local Nusselt number for each station Nu. 

(a) Enter figure 18(d) with My and (3 and determine 

Nu/ v^Ee. 

(b) The local Nusselt number for each station is then found 

by multiplying this value by the square root of the 
Reynolds numbers for the respective stations. 

16. Local heat— transfer coefficient h. 

(a) Enter figure 18(e) with T B ‘— 460 and find the thermal 

conductivity of air at the surface of the cone k_ ’ . 

Nu k ‘ 

(b) Then, h is equal to — B — for each station. 

17 . Local rate of heat transfer q. 


18. Average heat— transfer coefficient h. The average coeffi- 
cient for that portion of the cone from the apex to any point along 
its surface for laminar flow is given by h = (4/3)h. 
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Figure i Design details of the 20° cone. 
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Figure 2.— Electrically heated 20° cone with power terminals , voltage- 
tap leads and thermocouple leads. 



NACA 


Figure 3.- Heated 20° cone installed in the test section of the Ames 
1_ by 3— foot supersonic vind tunnel No. 1. 
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Figure 4. - Simplified wiring diagram for electrically heated cone. 
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Cone length ; s, in. 

(a) H 0 - 3.1 tb/sq in. abs. 

Figure 5. — Surface - temperature distributions for various 
nominal surface temperatures on the heated 20° 
cone with a laminar boundary layer. 
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Cone length , 5, in 
(b) H 0 s 5.7 !b/sq In.obs . 


Figure 5 - Continued . 


CONFIDENTIAL 




CONFIDENTIAL 


NACA EM No. A8L28 



0 2 4 6 8 

Cone length , s,/n. 

(c) H 0 - 8.8 lb/sq in. abs. 

Figure 5 -Continued. 
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Cone length, s , in. 
(d) Ho - 14.8 lb/sq in.abs. 

Figure 5r Continued. 
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Cone length , s, in. 

(e) H 0 S 2L0 lb/sq in. abs. 

Figure 5. - Concluded. 
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Local heat -transfer rate, q , Btu/hrft 



(a) H 0 = 3.1 !b/sq in. abs. 


Figure 6. — The variation of the local rate of heat transfer with surface - temperature 
parameter for successive positions along the length of the 20° cone. 
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Figure 6. — Continued. 
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(c) H 0 =8.8 lb/sq in. abs. 

Figure 6. — Continued. 
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(d) H 0 =14.8 /b/sq in. abs. 


Figure 6. — Continued. 
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Local heat -transfer rate, q, Btu/hrff 
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(e) H 0 -21.0 !b/sq in. abs. 


Figure 6 . — Concluded. 
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Figure 7.-4 comparison of the results obtained with and without 
the assumption of a negligible surface -temperature gradient. 
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Figure 8 . — Diagram of the method of correcting 
the laminar boundary - layer thickness for the 
effect of the unheated nose of the 20 cone . 
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Figure 9. —A comparison of experimental laminar heat- 
transfer-coefficient distribution on a 20° cone and the 
theoretical distributions for both the heated and unheated 
nose conditions. 
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(a) <3-14. 

Figure 10. — The variation of local Nusselt number with length 
Reynolds number for a laminar boundary layer on the 
heated 20 9 cone. 
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Figure !0 —Continued. 


CONFIDENTIAL NACA EM No. A8L28 


CONFIDENTIAL 



(c) /3 s 1.8 . 


Figure !0. — Continued. 


VJl 

uo 


NACA EM No. A8L28 CONFIDENTIAL 


CONFIDENTIAL 


VJI 

-p- 



(d) /3 -- 2.0. 


Figure 10 . — Concluded. 
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Figure II. — A comparison of the theoretical laminar heot- 
transfer— coefficient distributions on a 20° cone with 
the experimental distributions. 
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Figure 12. — Comparison of the experimental local heat- 
transfer-rate distribution with that calculated by 
the theory for laminar boundary layers. 
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Figure 13. ~ The variation of the length Reynolds number 
at transition with surface-temperature parameter on 
a heated 20° cone. 
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Figure 14.— Surface — temperature distributions for various 
nominal surface temperatures on the heated 20° 
cone with a turbulent boundary layer. 
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(b) H 0 *14.8 Ib/sq in . abs. 

Figure !4 .—Concluded. 
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Figure 15- Distribution of local rate of heat transfer 
on a 20° cone with a turbulent boundary layer. 
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Figure 1 5 .—Concluded . 
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a turbulent boundary layer. 
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Figure 17. —A comparison of experimental data with 
Eber's equation and the Hantzsche and Wendt 
equation tor a 20° cone with a turbulent boundary 

layer. 
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(a) Relation between surface Mach number and cone half -angle for various free- stream 
Mach numbers (reference II). 

Figure 18 . -Design charts. 
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(b) The variation of the ratio of recovery to total 
temperature with Mach number. 


Figure ! 8. — Continued. 
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(c) The variation of absolute viscosity of air with 

temperature. 

Figure 18. — Continued 
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(d) The variation of local Nusselt number with surface - 
temperature parameter , Mach number , and Reynolds 
number for cones with laminar boundary layers. 



Figure 18. — Continued. 
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(e) The variation of thermal conductivity of air with 
temperature. 

Figure 18. — Concluded . 
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